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Abstract. Classical physics is about real objects, like apples falling from trees, whose motion is 
governed by Newtonian laws. In Standard quantum mechanics only the wave function or the results 
of measurements exist, and to answer the question of how the classical world can be part of the 
quantum world is a rather formidable task. However, this is not the case for Bohmian mechanics, 
which, like classical mechanics, is a theory about real objects. In Bohmian terms, the problem of 
the classical limit becomes very simple: when do the Bohmian trajectòries look Newtonian? 

1 Introduction 

The classical world, say the world of objects of familiar experience that obey Newtonian laws, 
seems far removed from the "wavy" world of quantum mechanics. In this paper we shall 
sketch what we believe are the bàsic steps to be taken in going from the quantum world to 
the classical world. 

1. The first step is the crucial one: As is well known, and as Bell has emphasized [p]], 
Standard quantum mechanics is not a precise microscopic theory because the division between 
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the microscopic and the macroscopic world, which is essential to the very formulation of 
that theory, is not made precise by the theory ||. In fact, the following conclusion seems 
inevitable: quantum mechanics does not contain the means for describing the classical world 
in any approximate sense and one needs to go beyond quantum mechanics in order to do 
so. There are two natural possibilities for amending ordinary quantum mechanics: either 
the wave function is not all there is, or Schròdinger's equation is wrong. In this paper we'll 
formulate the problem of the classical limit within the framework of Bohmian mechanics, a 
theory which follows the first path and in which the observer doesn't play any crucial role. 
It is a theory about reality, not about the result of measurements. A very short review of 
Bohmian mechanics is given in section |2| and the the relevant part of this theory, related to 
the classical limit, is discussed in section f|. 

2. To get a handle on a problem, one should first simplify it as much as possible. The 
complex motion of a macroscopic body can be drastically simplified by making some rather 
Standard approximations and reducing the problem to that of a "particle" moving in an 
external potential. This is what we shall do in section |3|. 

3. Good textbooks on quantum mechanics contain enlightening ideas. One of these ideas is 
the so called Ehrenfest theorem that we shall use in section [5] in order to obtain a necessary 
condition for the classicality for wave packets. 

4. The structure of Bohmian mechanics contains the means for extending the condition for 
classicality to more general wave functions, namely wave functions which locally look like 
plane waves, as we shall see in section ^. In section [7] we shall then show how the problem of 
classical limit for general wave functions can be reduced to that for local plane waves. 

5. Simplicity is good, but it has its limitations: the reduction of the motion of the center of 
mass to a one body problem doesn't explain the robustness and stability of classical behavior. 
This however can be explained by making the model a little more realistic, say by including 
in an effective way the external as well as the infernal environment. We shall briefly touch 
this point in section || 

6. This step is the crucial one from a mathematical point of view. In section |9| we shall put 
forward a mathematical conjecture on the emergence of classical behavior. Unfortunately we 
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cannot provide any rigorous mathematical justification for it. Mathematical work on it would 
be valuable since this conjecture goes beyond the Standard mathematical work of semiclassical 
analysis (see, e.g., |L3|], |l·|) or, in more modern terms, microlocal analysis (see, e.g., JL2]). 



7. This is the last step in what we believe is the main structure of the classical limit: 

{iP,X)^(P,X), 

where on the two sides of the arrow are represented the complete state description of Bohmian 
mechanics, in terms of wave function and position, and of classical mechanics, in terms of 
momentum and position. 



2 Bohmian Mechanics 

In nonrelativistic Bohmian mechanics the world is described by partides which follow tra- 
jectòries determined by a law of motion. The evolution of the positions of these partides is 
guided by the wave function which itself evolves according to Schròdinger's equation. In other 
words, in Bohmian mechanics the complete description of the state of an iV-particle system 
is the pair Q), where ^ = ^(q) = ^(çi, ■ ■ ■ , çn) and Q = (Qx, ...,Qn) are respectively the 
wave function and the actual configuration of the system, with Qk denoting the position of 
the k-th particle in ordinary three-dimensional space. 

For non relativistic spinless partides the state Q) evolves according to the equations 

d(h h T v gfc g(Q) m 

dt m k *(Q) 1 1 

d^> N fi 2 „ 

Equations (HD and @ form a complete specification of the theory. Agreement between 
Bohmian mechanics and quantum mechanics regarding the results of any experiment is guar- 
anteed by what has been called H] the quantum equílíbríum hypothesis: when a system has a 
wave function if>, its configuration Q is random with probability distribution 

p(q) = mq)\ 2 . (3) 

While the meaning and justification of this hypothesis is a delicate matter, which has been 
discussed at length elsewhere H (see also Q and @), we wish to underline here an important 
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property of (0): if the probability density for the configuration satisfies p(q,t ) = \ip(q,t )\ 2 
at so me time to, then the density at any time t to which this is carried by the motion (0) is 
also given by p(q,t) = í)| 2 . This is an extremely important property of any Bohmian 
system, expressing a compatibility between the two equations of motion (P and (Q) defming 
the dynamics, which we call the equivariance of |"0| 2 - 



3 Motion in an External Potential 

Our goal is to study the classical behavior of a macroscopic body composed of N partides 
with N ^> 1 (one may think of an apple falling from a tree or a planet moving around the 
sun). It is rather clear that one expects classical behavior only for appropriate macroscopic 
functions of the particle configuration (Qi, ...,Qn). The relevant macroscopic variable, whose 
classical behavior we wish to investigate here, is the center of mass of the body 

x = Ei rriiQi 



m 

where mi, . . . , tun are the masses of the partides composing the body and m = J2i m i is the 
total mass of the body 

We shall assume that the partides interact through internal forces as well as being subjected 
to an external potential, so that the potential energy in (^) is of the form 

U(q)=J2U(q hqj )+J2Vi(qi). 

i<j i 

Let y = (yi, . . . ,yN-i) be a suitable set of coordinatesQ relative to the center of mass x = 
J2i rriiqi/m. Then under the change of variables q = (x, y) Schròdinger's equation (0) assumes 
the form 

ih^- = (h x + H y + H {x ' y) ) V (4) 

where 

H x = t^V, 2 + V(x) , V{x) = ^V t (x) } 
2m 

i 

H y is the free Hamiltonian associated with the relative coordinates y and the operator H^ x,v ^ 
describes the interaction between the center of mass and the relative coordinates. If V; are 



1 For sake of concreteness one may think, e.g., of the so called Jacobi coordinates. 
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slowly varying on the size of the body, H( x > y > can be treated as a small perturbation, and, in 
first approximation, neglected. Thus, if \I/ = ip(x)(f)(y) at some time, the time evolution of 
the center of mass decouples from that of the relative coordinates and we end up with a very 
simple one particle problem: the wave function ip of the center of mass evolves according to 
one-particle Schròdinger's equation 

ih JÏ = l~ V ^ + ( 5 ) 
at 2m 

and its position X evolves according to 

dx h v^(x) ... 

From now on, whenever no ambiguity will arise, we shall treat the center of mass as a "particle" 
and we shall refer to X and if> as the position and the wave function of such a particle. 



4 The Classical Limit in Bohmian Mechanics 

In order to investigate the conditions under which X evolves classically it is useful to write 
the wave function if) = ifj(x) in the polar form 

V>(a;) = R(x)eï s{x \ (7) 

From Schròdinger's equation ([5]) one obtains, following Bohm , the continuity equation for 



m 

~dt 



2 



+ div 



V a 5 



m 



R 



0. 



and the modified Hamilton- Jacobi equation for S 



dt 2m 2m R w 

Note that equation (|S|) is the usual classical Hamilton- Jacobi equation with an additional 
term 

h 2 V X 2 R . , 

v ° s -sí-s-· (10) 

called the quantum potential. Since is the right hand side of @, one then sees that 
the (size of the) quantum potential provides a rough measure of the deviation of Bohmian 
evolution from its classical approximation. 
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Analogously, consider the modified Newton equation associated with (0), and obtained by 
differentiating both sides of equation @ with respect to time, 

m ^ = F + i7 Q' ( n ) 

where F = — V X V(X) and Fq = — V x Vq(X) are respectively the classical force and the 
" quantum" force. Equation (|ÏI]) shows that all the devíations from classicality are embodied 
in the quantum force Fq. 

Thus, the formulation of the classical limit in Bohmian mechanics turns out to be rather 
simple: classical behavior emerges whenever the particle trajectory X = X(t), satisfying (|TT|), 
approximately satisfies the classical Newton equation, i.e., 

dt 2 v ; 

The problem is to determine the physical conditions ensuring (0). Usually, physicists consider 
classical behavior as ensured by the limit h —>■ 0, meaning by this 

h < A , (13) 

where Aq is some characteristic action of the corresponding classical motion (see, e.g., 
|^3[1,|Ï5|,P1)). Condition ( |Ï3"D is often regarded as equivalent to another Standard condition 



of classicality which involves the length scales of the motion (see, e.g., fl p): if the de Broglie 
wave length À is small with respect to the characteristic dimension L determined by the scale 
of variation of the potential V, the behavior of the system should be close to the classical 
behavior in the same potential V. This is very reminiscent of how geometrical òptics can be 
deduced from wave òptics. We regard this condition, i.e., 

A < L, (14) 

as the most natural condition of classicality since it relates in a completely transparent way 
a property of the state, namely its de Broglie wave length À, and a property of the dynamics, 
namely the scale of variation of the potential L. In the remainder of this paper we shall argue 
that (|H] ) is indeed a necessary and sufficient condition for (|T^). 
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5 Wave Packets 



To explain the physical content of ( H ) and its implications we shall consider first the case 
for which the wave function has a well-defined de Broglie wave length: we shall assume that 
^ is a wave packet with diameter a, with mean wave vector k and associated wave length 
A = 2vr/|fc|. 

As we shall see, the analysis of this situation will allow us to find a precise characterization 
of the scale L of variation of the potential. Our analysis will be rather Standard — it is basically 
the Ehrenf est's Theorem — and can be found in good textbooks (see, e.g., jÏQfl). We reproduce 
it here both for the sake of completeness and because we believe that it attains, within 
the Bohmian framework, a deeper and much more general significance than within Standard 
formulations of quantum mechanics. 

From the equivariance of (||D we have that the mean particle position at time t is given by 

(X) = / x\^ t (x)\ 2 dx. 



From (^|) it follows that 

m ^L(X) = -J V x V{x)\^ t {x)\Hx. 
By expanding F(x) = — V x V{x) in Taylor series around (X) one obtains 



dxjdxk 



(15) 



where 



Ajfi — (XjX k ) — (Xj)(X k ) 

is of order cr 2 , where a is the diameter of the packet. Therefore, the mean particle position 
should satisfy the classical Newton equation whenever 



d 3 V 



dxidxjdxk 



dV 



dxi 



í.e., 



a < 



V 



(16) 



;i7) 



where V and V" denote respectively suitable estimates of the first and third derivatives (e.g., 
by taking a sup over the partial derivatives). 
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The minimum value of the diameter of the packet a is of order À. Hence (O) becomes 

A < 



V 



v 

This last equation gives a necessary condition for the classicality of the particle motion and, by 
comparing it with (PS), a precise definition of the notion of scale of variation of the potential, 
namely, 



L(V) 



V 



V'- 



(19) 



In the following we shall argue that (0), with L given by (|Ï9"D, is indeed also sufficient 
for classical behavior of Bohmian trajectòries. For wave packets this follows easily from 
the equivariance of \ip\ 2 : over the lapse of time for which the spreading of the packet can 
be neglected, the overwhelming majorityQ of trajectòries X = X(t) will stick around their 



mean value (X) and follow its classical time evolution. Thus we expect ( fï2|) to hold for the 
overwhelming majority of trajectòries. 



6 Local Plane Waves 

Suppose now that ip is not a packet but a wave function that locally looks like a packet. By 
this we mean, referring to the polar representation (^), that the amplitude R(x) and the local 
wave vector 

k = k(x) = V x S(x)/h (20) 

are slowly varying over distances of order À(x) = h/\'V x S(x)\, the local de Broglie wave length. 
We may call such a ip a "local plane wave" . 

At any given time the local plane wave can be thought as composed of a sum of wave 
packets: Consider a partition of physical space into a union of disjoint sets Aj chosen in such 
a way that the local wave vector k(x) doesn't vary appreciably inside of each of them and 
denote by ki the almost constant value k(x) for ï 6 Aj. Let xAí be the characteristic function 
of the set Aj (xa z ( x ) = 1 if £ G Aj and otherwise). Since J2íXAí = 1, we have 

^( x ) = Ç xaA x )^( x ) = Ç ^i( x )- ( 21 ) 

i i 

2 With respect to the equivariant measure \ip\ 2 - 
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Note that this decomposition is somewhat arbitrary: provided that k(x) is almost constant 
in Aj, the extent of these sets can be of the order of many wave lengths down to a minimal 
size a-, ~ (A^l 1 / 3 of the same order of À ? n 

At any time, the position X of the particle will be in the support of one of the packets 
forming the decomposition (21), say in the support of fy. If the condition flT7p holds for <7j, 
we may then proceed as in the previous section: the minimal size of the packet ipi can be 
taken of order <jj = X(x) and the condition of classicality is again (18) for A = X(x). 

Note that this straightforward reduction of the classical limit for local plane waves to that 
for wave packets is possible only within Bohmian mechanics: since the particle has at any 
time a well-defmed position X and the different components of the local plane wave ( pTf ) don't 
interfere, we may "collapse" ip to the wave packet ipi relevant to the dynamics of X . 



7 General Wave Functions 

We wish now to investigate the physical content of (fH| ) and its implications for a general 
wave function. The first issue to address is what notion of wave length should be appropriate 
for this case. A rough estimate of À could be given in terms of mean kinetic energy associated 
with ip, 

h 2 

E Mn (iP) = (iP,-—V x 2 iP), (22) 

with associated wave length 

A = À(V0 = 1=^=== • (23) 

Suppose now that (|ï"8|), with À given by (p3|), is satisfied. We claim that in this case the 
Schròdinger evolution should "quickly" produce a local plane wave, that can be effectively 
regarded as built of pieces that are wave packets satisfying ( |Ï8| ) for À = À(x) and hence 
themselves evolving classically as we have seen in the previous section. 

In fact, if À <C L the kinetic energy dominates the potential energy and the free Schròdinger 
evolution provides a rough approximation of the dynamics up to the time needed for the 
potential to affect the evolution significantly. During this time, the Schròdinger evolution 
produces a spatial separation of the different wave vectors contained in tp, more or less as 

3 The use of the characteristic function may introduce an undesirable lack of smoothness, but this can be 
easily taken care by replacing the \Aí with functions Oi forming a smooth partition of unity. 
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Newton's prism separates white light into the different colors of the rainbow. In other words, 
the formation of a local plane wave originates in the " dispersive" character of free Schròdinger 
evolution. 

So, in order to gain some appreciation of this phenomenon consider the free Schròdinger 
evolution 

1 r it 



, x hk 



Mx) = j2^J e ' ""W)**, ( 24 ) 

where $ is the Fourier transform of the initial wave function ip. The stationary phase method 
yields straightforwardly the long time asymptotics of ipt, 

/ 2, i yY} / \ m ^2 A 777 CC 

ipt(x) ~ í t - ) Q l ^~il){k) , where = — — , (25) 

which is indeed a local plane wave with local wave vector k = mx/(ht). 

We said above that the local plane wave is "quickly" produced. But how quickly? In order 
to estimate such a time, consider the simple example of an initial wave function ip composed of 
two overlapping wave packets with the same position spread Ax and with opposite momenta 
p and —p. The time r of formation of a local plane should be of the order of the time for 
separation of the packets, which is basically the time needed to cover a space equal to Ax. 
From AxAp ~ H and Ap ~pwe obtain 

Ax h h 

T „ (26) 

p/m p z /m {E} 

where (E) is the mean kinetic energy of the particle. It is reasonable to suggest that (|26|). 
with (E) given by (|2~2"D , could give a very rough estimate of the time of formation of a local 
plane wave for a general wave function ip. Note that the time needed for the potential to 
produce significant effects on the evolution is of order 

T = — , where v = — -. (27) 

Thus, if À L we have that t ^ T, which means that the local plane wave gets formed on 
a time scale much shorter than the time scale over which the potential affects the dynamics. 

We arrive in this way at a sharp (or, at least, sharper than usually encountered) mathe- 
matical formulation of the classical limit for a general wave function ip. First of all, consider 
the dimensionless parameter 

e - A W Í281 
e -L(V)- (28) 
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Secondly, consider the Bohm motion X on the "macroscopic" length and time scales defined 
by ip and V. By this we mean X' = X'(t'), where 



X' = X/L and t' = t/T (29) 

with T given by (0). Finally, consider Fq/tu, the "quantum" contribution to the total 
acceleration in ([ïï]) , on the macroscopic scales (p9j), namely 

D = ^F Q (X'L,t'T) (30) 

Then the Bohm motion on the macroscopic length and time scales will be approximately 
classical, with deviation from classicality D tending to as e — ► 0. 



We'd like to point out that the use of macroscopic coordinates ( |29| ) for the formulation 
of the classical limit is rather natural from a physical point of view. First of all, the scales 
L and T are the fundamental units of measure for the motion: L is the scale on which the 
potential varies and T provides an estimate of the time necessary for the particle to see its 
effects. More importantly, in the limit e — > the nonclassical behavior — occurring during the 
time r of formation of the local plane wave — disappears, since, as we have argued above, in 
this limit t <C T. In other words, on the macroscopic scales on which we expect classical 
behavior the local plane wave has been formed. 



8 Limitations of the Model: Interference and the Role 
of the Environment 

Before commenting on the mathematics of the limit e — ■> we should stress a physical caveat. 
For motion in unbounded space, the expanding character of the Schròdinger evolution makes 
the set of local plane waves an "attractor" for the dynamics — so that the local plane wave 
form is in this sense "typical" . However, for motion in a bounded region (with wave functions 
which are superpositions of bound states) the "typical" wave function is composed by a sum 
of local plane waves, this being due to interference between the waves reflected by the "edges" 
of the confining potential. Consider for example an infmite potential well of size L in one 
dimension and initial wave function ip, well localized in the center of the well which is the 
superposition of two packets with opposite momenta p and —p. Suppose that X(ip) <C L. 
Then the two packets move classically and at a certain time, say t r , are reflected from the 



11 



walls of the potential. At the time t c = 2t r , they interfere in the middle of the well. t c is the 
u first càustic time," the time at which the classical action S c \(x,t) becomes multivalued. In 
general, we should not expect classical behavior for times larger than the first càustic time t c . 

What is going on? The emergence of classical behavior should be robust and stable, which 
would not be the case if it were restricted to times smaller than t c . However, if one remembers 
that the model we are investigating is a strong idealization, the problem evaporates. We 
are in fact dealing with the one-body problem defined by (§) and (|6|), an approximation 
to the complete dynamics defined by (§) in which the term H^ x ' v \ describing the interaction 
between the center of mass and the relative coordinates, is neglected. Note than even (^) is an 
idealization since it does not include the unavoidable interaction of the body with its external 
environment: in a more realistic model H^ x ' y ' would take into account both the infernal and 
external environment of the center of mass (with y now including both the relative coordinates 
and the degrees of freedom of the external environment). These interactions — even for very 
small interaction energy — should produce entanglement between the center of mass x of the 
system and the other degrees of freedom y, so that their effective role is that of "measuring" 
the position X and suppressing superpositions of spatially separated wave functions. (Taking 
these interactions into account is what people nowadays call decoherence, see, e.g., |7j and 
the references therein). Referring to the above example, the effect of the environment should 
be to select (as relevant to the dynamics of X, see f| and [|]) one of the two packets on a 
time scale much shorter than the first càustic time t c . 

9 Towards a Mathematical Conjecture 

The mathematical content of sections |7] and |8| is summarized by the following (not yet sharply 
formulated) conjecture: 

Conjecture. Let e be the dimensionless parameter defined by (ffíBQ and D be the quantity given 
by ffify- Then there are environmental interactions such that D ~> as e —>■ 0, uniformly in 
ip and V. 

Concerning this conjecture, we'd like to make here just a few remarks. 

1. U D ~> 0" means convergence to in a "suitable" probabilistic sense since D is a random 
variable. D is a function of X, and X is random with probability distribution given by 
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\ip\ 2 . To require almost sure convergence is probably too strong a demand. Convergence in 
probability, or L 2 convergence, would seem more appropriate. Moreover, À is defined in (p3|) 
in terms of the average kinetic energy. This average could be large even when there is a 
significant probability for a very small kinetic energy. Thus it is probably necessary to regard 
À as random (with randomness inherited from the kinetic energy) and to understand e — > 
also in a probabilistic sense. 

2. Uniformity of the limit in ip and V could be expressed as follows: let (V n ,íp n ) be any 
sequence for which e n = ^ — > 0, with \ n = X(ip n ) given by (|23|) , and L n = L(V n ). Then 
D ~» as n — > +oo. Understanding D ~> as convergence in probability, we could also 
express uniformity in the following way: for any r] > and for any 5 > 0, there exists an 
€q > such that P(D > 5) is smaller than rj whenever e < 6q- Here P is the probability 
measure defined by l^] 2 , i.e., P(dq) = \^/(q)\ 2 dq, which includes randomness arising from the 
environment. 

3. For quadratic potentials (including free motion and motion in a uniform force field) L = oo 
so that e = 0. In this case the conjecture should be modified as follows: let L Q be any length 
scale and T Q the corresponding time scale T Q = Then for the Bohm motion on the scales 
given by L a and T a , D ~> uniformly in ib and L D whenever e = — > 0. 

4. There is an enormous amount of mathematical work, called semiclassical analysis or, in 
more modern terms, microlocal analysis, in which the limit H — > of Schròdinger evolutions 
is rigorously studied. It should be stressed that the limit e — > is much more general than the 
limit h — > 0. In fact e — X/L — h/mvL. So keeping L and the momentum mv fixed, the limit 
h — > implies e — > 0. But there are many ways in which e could go to zero. The classical 
limit, as expressed by the above conjecture, is (at the very least) a two-parameters limit, 
involving À and L, and h —>■ is just a very special case. Moreover, these two parameters 
themselves live on infinite dimensional spaces since À = A (■?/>), with ip varying in the Hilbert 
space of the system's wave functions, and L = L(V), with V varying in the class of admissible 
one particle potentials (that is, potentials leading to a self-adjoint Hamiltonian) . 

5. Exactly for the reason expressed in the previous remark, the conjecture is really very hard 
to prové: it require a lot of uniformity both in the wave function ip and in the potential V. 
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Just to have an idea of the difficulties, one may think of the analogous problem in statistical 
mechanics, namely the problem of studying the deviations from thermodynamic behavior of 
a large but finite system about which not so much is known. 

6. While the conjecture is difhcult to prové, it is still not completely satisfactory from a 
physical point of view. The conjecture states only that D depends on e in such a way that 
D ~» as e — > 0, uniformly in ip and V . A physically more relevant result would be to 
estimate how rapidly D is tending to (e.g., like e, or e 2 or whatever). Note that only this 
last kind of result can be of practical value: given V and ip, it provides an estimate for the 
deviation from classicality, while any other results do not quite do this. 

7. With the conjecture, and even with the refmement proposed in the previous remark, 
there is a further difficulty to consider: even if H^ x ' v ' is treated as a small perturbation in 
the suggestion of section [8] might not be too realistic. In fact, the autonomous Schròdinger 
evolution, even of a very narrow wave function ip = ip(x), could be destroyed in very short 
times. This is a serious difficulty; one resolution might be found in the notion of condítional 
wave function of the x-system, ip(x) = *&(x, Y), where Y is the actual configuration of the 
environment (this notion has been introduced and analyzed in 0). We regard the extension 
of the conjecture to this more realistic framework as the most interesting open problem on 
the classical limit, which we leave for future work. 

10 The Classical Limit in a Nutshell 

The key ingredient in our analysis of the emergence of the classical world, is that as soon as 
the local plane wave has formed, each configuration X is attached to a guiding wave packet 
with a definite wave vector k(x,t) that locally determines the particle dynamics according to 
the local de Broglie relation 

p(x, t) = hk(x, t), 

which, for À <C L, evolves according to classical laws. This means that the classical limit can 
be symbolically expressed as 

where (ip, X) is the complete quantum state description in terms of wave function and position, 
while (P, X) is the complete classical state description in terms of momentum and position. 
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All the relevant macroscopic information contained in the pair (ip, X) is, in the classical limit, 
embodied in the pair (P,X) — the only robust, stable quantity. In other words, as far as the 
macroscopic dynamics of X is concerned, only the information carried by P is relevant. 
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